EXPLICIT RATES OF APPROXIMATION 
IN THE CLT FOR QUADRATIC FORMS 

FRIEDRICH GOTZEi AND ANDREI YU. ZAITSEV^ ^ 

Abstract. Let X, Xi, X2, ■ ■ ■ be i.i.d. R'^-valued real random vectors. Assume that 
EX = , cov X = C, E llXlp = and that X is not concentrated in a proper subspace 
of M''. Let G be a mean zero Gaussian random vector with the same covariance operator 
as that of X. We study the distributions of non-degenerate quadratic forms (Q[S'Ar] of the 
normahzed sums Sn — N^^^^ {Xi + • • • + Xn) and show that, without any additional 
conditions, 



A dof 
L^N — sup 



[^at] <x}-P{Q[G] <a;} =0{N-^), 



provided that c? > 5 and the fourth moment of X exists. Furthermore, we provide 
explicit bounds of order 0(iV~^) for An for the rate of approximation by short asymp- 
totic expansions and for the concentration functions of the random variables Q[S'jv + a], 
a € R*^. The order of the bound is optimal. It extends previous results of Bentkus and 
Gotze (1997a) (for d > 9) to the case d > 5, which is the smallest possible dimension 
for such a bound. Moreover, we show that, in the finite dimensional case and for iso- 
metric Q, the implied constant in 0{N-^) has the form Cd.a'^ (det C)-^/^ E IjC-^/^ X]]'' 
with some Cd depending on d only. This answers a long standing question about optimal 
rates in the central limit theorem for quadratic forms starting with a seminal paper by 
C.-G. Esseen (1945). 



1. Introduction 

Let M.'^ be the (i-dimensional space of real vectors x = {xi, . . . , Xd) with scalar product 
{x,y) = xiUi + ■■■ + XdUd and norm ||x|| = (a;, x)^/^. We also denote by a real 
separable Hilbert space consisting of all real sequences x = {xi,X2, ■ ■ ■) such that = 
xf + xl + ■ ■ ■ < 00. 

Let X, Xi, X2, ■ . . be a sequence of i.i.d. M'^-valued random vectors. Assume that 
EX = and = E \\Xf < 00. Let G be a mean zero Gaussian random vector such 
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that its covariance operator C = cov G : M*^ — M*^ is equal to covX. It is well-known that 
the distributions C{Siy) of sums 

S^''£N-'/'{Xr + --- + XN) (1.1) 

converge weakly to C{G). 

Let Q : — 7- M°' be a linear symmetric bounded operator and let Q[x] = {Qx,x) be 
the corresponding quadratic form. We shall say that Q is non-degenerate if ker Q = { } . 

Denote, for q > 0, 

(3/=B\\xr, (3 = (3,. 
Introduce the distribution functions 

F{x) = P{Q[S^] < x], H{x) '=^'P{Q[G] < x}. (1.2) 

Write 

An = sup\F{x)-H{x)\. (1.3) 

Theorem 1.1. Assume that Q and C are non- degenerate and that d > 5 or d = oo. 
Then 

An<c{Q,C)(3/N. 
The constant c(Q, C) in this bound depends on Q and C only. 

Theorem 1.2. Let the conditions of Theorem \l.l\ be satisfied and let 5 < d < oo. Assume 
that the operator Q is isometric. Then 

Atv < Q cr"^ (det C)-^/^ E \\C-'/^ XW^N. 

The constant Cd in this bound depends on d only. 

Theorems ll.ll and ll.2l are simple consequences of the main result of this paper, Theo- 
rem [2?3] (see also Theorem 12. ip . Theorem 11.11 was proved in Gotze and Zaitsev (2008). It 
confirms a conjecture of Bentkus and Gotze (1997a) (below BG (1997a)). It generalizes 
to the case d > b the corresponding result of BG (1997a). In their Theorem 1.1, it was 
assumed that (i > 9, while our Theorem 11.11 is proved for d > 5. Theorem 11.21 yields an 
explicit bound in terms of the distribution C{X). 

The distribution function of HS'^rp (for bounded X with values in W^) may have jumps 
of order (9(A^~^), for all 1 < ci < cxd. See, e.g., BG (1996, p. 468). Therefore, the bounds 
of Theorems 11.11 and 11.21 are optimal with respect to the order in N . 

Theorems \\.\\ \m\ and the method of their proof are closely related to the lattice point 
problem in number theory. Suppose that d < oo and that (Qx, x) > 0, for x ^ 0. Let 
vol Er be the volume of the ellipsoid 

Er = {x eM.'^ : Q[x] <r^}, r > 0. 

Write volz Er for the number of points in ErdZ'^, where Z"' C M'^ is the standard lattice 
of points with integer coordinates. 
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The following result due to Gotze (2004) is related to Theorems 11.11 and 11.21 (see also 
BG (1995a, 1997b)). 

Theorem 1.3. For all dimensions d > 5, 



sup 



volz {Ej. + a) — vol Er 

VolEr 



0{r'^), for r > 1, 



where the constant in 0{r ^) depends on the dimension d and on the lengths of axes of 
the ellipsoid Ei only. 

Theorem 11.31 solves the lattice point problem for d > 5. It improves the classi- 
cal estimate C)(r~^'^/*^'^+^^) due to Landau (1915), just as Theorem 11.11 improves the 
bound C(iV-'^/('^+^)) by Esseen (1945) in the CLT for ellipsoids with axes parallel to 
coordinate axes. A related result for indefinite forms may be found in Gotze and Mar- 
gulis (2010). 

Work on the estimation of the rate of approximation under the conditions of Theo- 
rem 11.11 for Hilbert spaces started in the second half of the last century. See Zalesskii, 
Sazonov and Ulyanov (1988) and Nagaev (1989) for optimal bounds of order 0{N~^^'^) 
(with respect to eigenvalues of C) assuming finiteness of the third moment. For a more de- 
tailed discussion see Yurinskii (1982), Bentkus, Gotze, Paulauskas and Rackauskas (1990), 
BG (1995b, 1996, 1997a) and Senatov (1997, 1998). 

Under some more restrictive moment and dimension conditions the estimate of order 
C(A^-^+"), with £ ; as (i t oo, was obtained by Gotze (1979). The proof in Gotze (1979) 
was based on a new symmetrization inequality for characteristic functions of quadratic 
forms. This inequality is related to Weyl's (1915/16) inequality for trigonometric sums. 
This inequality and its extensions (see Lemma 16. ip play a crucial role in the proofs of 
bounds in the CLT for ellipsoids and hyperboloids in finite and infinite dimensional cases. 
Under some additional smoothness assumptions, error bounds 0{N~^) (and, moreover, 
Edgeworth type expansions) were obtained in Gotze (1979), Bentkus (1984), Bentkus, 
Gotze and Zitikis (1993). BG (1995b, 1996, 1997a) established the bound of order O^N'^) 
without smoothness-type conditions. Similar bounds for the rate of infinitely divisible 
approximations were obtained by Bentkus, Gotze and Zaitsev (1997). Among recent 
publications, we should mention the papers of Nagaev and Chebotarev (1999), (2005) 
{d > 13, providing a more precise dependence of constants on the eigenvalues of C) and 
Bogatyrev, Gotze and Ulyanov (2006) (non- uniform bounds for d > 12), see also Gotze 
and Ulyanov (2000). The proofs of bounds of order 0{N~^) are based on discretization 
(i.e., a reduction to lattice valued random vectors) and the symmetrization techniques 
mentioned above. 

Assuming the matrices Q and C to be diagonal, and the independence of the first 
five coordinates of X, Bentkus and Gotze (1996) have already reduced the dimension 
requirement for the bound 0{N~^) to d > 5. The independence assumption in BG (1996) 
allowed to apply an adaption of the Hardy-Littlewood circle method. For the general case 
described in Theorem II. H one needs to develop new techniques. Some yet unpublished 
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results of Gotze (1994) provide the rate 0{N^^) for sums of two independent arbitrary 
quadratic forms (each of rank d > 3). Gotze and Ulyanov (2003) obtained bounds of 
order 0{N~^) for some elhpsoids in M"^ with c? > 5 in the case of lattice distributions 
of X. 

The optimal possible dimension condition for this rate is just d > 5, due to the lower 
bounds of order 0{N~^ log N) for dimension d = 4 in the corresponding lattice point 
problem. The question about precise convergence rates in dimensions 2 < d < 4 still 
remains completely open (even in the simplest case where Q is the identity operator 1^, and 
for random vectors with independent Rademacher coordinates). It should be mentioned 
that, in the case d = 2, a precise convergence rate would imply a solution of the famous 
circle problem. Known lower bounds in the circle problem correspond to the bound of 
order 0{N~^/'^ log'^iV), 6 > 0, for A^. Hardy (1916) conjectured that up to logarithmic 
factors this is the optimal order. 

Now we describe the most important elements of the proof. We have to mention that 
a big part of the proof repeats the arguments of BG (1997a), see BG (1997a) for the 
description and application of symmetrization inequality and discretization procedure. In 
our proof we do not use the multiplicative inequalities of BG (1997a). Here we replace 
those techniques by arguments from the geometry of numbers, developed in Gotze (2004), 
combined with effective equidistribution results by Gotze and Margulis (2010) for suit- 
able actions of unipotent subgroups of SL(2,]R), see Lemma [8l2l These new techniques 
(compared to previous) results are mainly concentrated in Sections [SHHl 

Using the Fourier inversion formula (see (14.31) and (14.41) ). we have to estimate some 
integrals of the absolute values of differences of characteristic functions of quadratic forms. 
In Section [6l we reduce the estimation of characteristic functions to the estimation of a 
theta-series (see Lemma 1^751 and inequality (16.271) ). To this end, we write the expecta- 
tion with respect to Rademacher random variables as a sum with binomial weights p{m) 
and p{fn). Then we estimate p{m) and p{fn) from above by discrete Gaussian exponen- 
tial weights Csq{m) and Csqijn), see f l6.15p . f l6.18p . f l6.20p and f l6.2ip . Together with the 
non-negativity of some characteristic functions (see (I6.19P and (I6.23P ). this allows us to 
apply then the Poisson summation formula from Lemma 16.41 This formula reduces the 
problem to an estimation of integrals of theta-series. Section [7] is devoted to some facts 
from Number Theory. We consider the lattices, their a-characteristics (which are defined 
in (17. lip and fl7.12p ) and Minkowski's successive minima. In Section |8l we reduce the es- 
timation of integrals of theta-series to some integrals of a-characteristics. An application 
of the crucial Lemma 18.21 decribed above, ends the proof. 



2. Results 

To formulate the results we need more notation repeating most part of the notation 
used in BG (1997a). Let cr1> al> . . . be the eigenvalues of C, counting their multiplic- 
ities. We have a'^ = al + al + ■ ■ ■ . 
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We shall identify the linear operators and corresponding matrices. By 1^; : R'^ — )■ 
we denote the identity operator and, simultaneously, the diagonal matrix with entries 1 
on the diagonal. By we denote the {d x d) matrix with zero entries. 

Throughout 5 = {ei, . . . , e^} C M"^ denotes a finite set of cardinality s. We shall write 
So instead of S if the system {ei, . . . , e^} is orthonormal. 

Let p > and 6 > 0. Following BG (1997a), we introduce a somewhat modified 
non-degeneracy condition for the distribution of a ci-dimensional vector Y: 

Ar{p,6,S,Y) : P{\\Y -e\\<6}>p, for all e G 5. (2.1) 

We shall refer to condition (12. ip as condition Af{p, 6, S, Y). We shall write 

Mj(p, 6, S, Y) = ATip, 6, S, Y) U M{p, 6,QS, Y). 

Just condition Mq{p,6,S,Y) was used in BG (1997a). Note that 

Af{p,6,S,Y)=Af,^{p,6,S,Y). (2.2) 

Introduce truncated random vectors 

= XI{ ||X|| <(TyiV}, X^ = XI{\\X\\ > a^/N}, (2.3) 



= XI{\\C-^/^ X\\ <VdN}, Xa = XI{\\C-^^^ X\\ > VdN}, (2.4) 
and their moments (for q > 0) 

At=^E\\xr. Ii;=^-^E\\X4\ (2.5) 

A?= ^EllC-'/^Jfar, n°= -^-^E||C-"^Xor. (2.6) 

Here and below denotes the indicator of an event A. Of course, definitions ( 12. 4 p 

and (12. 6p have sense ii d < oo and the covariance operator C is non-degenerate. 
Clearly, we have 

X'' + X^ = X" + Xa=X, IIX^II ||X<>|| = ||X°|| llXnII = 0. (2.7) 

Generally speaking, X^ and X* are different truncated vectors. In BG (1997a) the i.i.d. 
copies of the vectors X* and X^ only were involved. Truncation (12. 4p was there applied 
to the vector X^. The use of X^ is more natural for the estimation of constants in the 
case d < oo. It is easy to see that 

(C-i/^X)" = (C-1/2X)° = C-i/2x°, (2.8) 

and 

(C-1/2X)^= (C-i/2^)^ = C-i/2^n. (2.9) 
Equalities (12. 8p and (12. 9p provides a possibility to apply auxiliary results obtained in 
BG (1997a) for truncated vectors X" and Xo to truncated vectors C^^/^ X° and C"^/^ X^. 
However, one should take into account that o"^, A4, H^, G, . . . have to be replaced by 
corresponding objects related to the vector C^^^"^ X (that is, by d, A4 , H°, C~^/^ G, . . . ). 
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In Sections m and we shall denote 

X' = X"-EX" + W, (2.10) 

where is a centered Gaussian random vector which is independent of all other random 
vectors and variables and is chosen so that covX' = covG. Such a vector W exists by 
Lemma 13.21 

By c, ci,C2,... we shall denote absolute positive constants. If a constant depends 
on, say, s, then we shall point out the dependence writing Cs or c(s). We denote by c 
universal constants which might be different in different places of the text. Furthermore, 
in the conditions of theorems and lemmas (see, e.g.. Theorem 12. H 12.21 and the proofs of 
Theorems 12. 3^ 12.71 and 12. 8p we write Cq for an arbitrary positive absolute constant, for 
example one may choose Cq = 1. We shall write A B, ii there exists an absolute 
constant c such that A < cB. Similarly, A B, if A < c{s)B. We shall also write 
A B if A -Cs B <t^s A. By \a\ we shall denote the integer part of a number a. 

Throughout we assume that all random vectors and variables are independent in ag- 
gregate, if the contrary is not clear from the context. By Xi,X2,... we shall denote 
independent copies of a random vector X. Similarly, Gi,G2, ■ ■ ■ are independent copies 
of G and so on. By C{X) we shall denote the distribution of X. Define the symmetrization 
X of a random vector X as a random vector with distribution C{X) = C{Xi — X2). 

Instead of normalized sums S^, it is sometimes more convenient to consider the sums 
Zj\f = Xi + ■ ■ ■ + Xtv. Then Sj^ = N'^^'^Z^. Similarly, by (resp. and Z'j^) 
we shall denote sums of N independent copies of X* (resp. X'^ and X'). For example, 
Z'j^ = X[ + ■ ■ ■ + X^. 

The expectation Ey with respect to a random vector Y we define as the conditional 
expectation 

Ey/(X,y,Z,...) = E(/(X,F,Z...)|X,Z,...) 

given all random vectors but Y. 

Throughout we write e{x} =^exp{ix}. By 

/oo 
e{tx}dF{x) (2.11) 
-00 

we denote the Fourier-Stieltjes transform of a function F of bounded variation or, in other 
words, the Fourier transform of the measure which has the distribution function F. 

Introduce the distribution functions 
Fa{x) = P{Q[SN-a]<x}, Ha{x) = P{Q[G-a]<x}, a G s G R. (2.12) 

Furthermore, define, for d = 00 and a G M^, the Edgeworth correction 



Ea{x) = E,{x;Q,C{X),C{G)) 
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as a function of bounded variation such that oo) = and its Fourier-Stieltjes trans- 
form is given by 

Ea{t)= ^^Ee{tQ[r]}(3(QX,F)(QX,X) + 2zt(QX,r)3), Y = G - a. (2.13) 

In finite dimensional spaces (for 1 < d < oo) we define the Edgeworth correction as 
follows (see Bhattacharya and Rao (1986)). Let denote the standard normal density 
in R'^. Then p{y) = (f){C-^/^y)/VdetC, y G is the density of G, and, for a G 
b = a/ZV a, we have 

E,(x) = eb(iVx) = -^x(^x), A^ = {ueR'': Q[u-a]<x}, (2.14) 
with the signed measure 

X{A)= [ Ep"'{y)X^dy, for the Borel sets A C (2.15) 

J A 

and where 

p'"{yW = p{y){3{C-\,u){C-'y,u) - {C-'y,uf) (2.16) 

denotes the third Frechet derivative of p in direction u. 

Notice that = if a = or if E (X, yf = 0, for all y G R'^. In particular, = if 
X is symmetric (that is, C{X) = £(— X)). 

We can write similar representations for E^{x) = B^(Xx), E^{x) = 0^(Xx) and 
E'^{x) = e'f,{Nx) just replacing X by X°, X^ and X' in fl27[3|) or f lfTSj) . 

For b G M'^, introduce the distribution functions 

^t{x) = P{Q[ZN-b]<x}=Fa{x/N), (2.17) 

and 

$fe(x) = PlQiViVG-fe] < a;} = Ha{x/N). (2.18) 
Define, for a G M'^, 6 = Vx a, 

Ai^) = sup I Faix) - Haix) - Ea{x) | = sup | ^ b{x) - ^,{x) - Q,{x) \ (2.19) 

(see ( I2l2|) . dm]) . ( 12X71) and ( l2l8|) to justify the last equahty in ( I2l9l) ). We write Aj^^^g 
and aS^^^^ replacing E„ by E^ and in f l219|l . 

The aim of this paper is to derive for is!"^ explicit bounds of order (9(X~^) without any 
additional smoothness type assumptions. Theorem 12.11 (which was proved in BG (1997a)) 
solved this problem in the case 13 < d < oo. 

In Theorems I2.1H2.8I we assume that the symmetric operator Q is isometric, that is, 
that is the identity operator 1^. This does not restrict generality (see Remark 1.7 in 
BG (1997a)). Indeed, any symmetric operator Q may be decomposed as Q = QiQoQi , 
where Qo is symmetric and isometric and Qi is symmetric bounded and non-negative, 
that is, (Qix,x) > 0, for all x G R'^. Thus, for any symmetric Q, we can apply all our 
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bounds replacing the random vector X by QiX, the Gaussian random vector G by QiG, 
the shift a by Qio, etc. In the case of concentration functions (see Theorems 12. 71 and 12.81) . 
we have Q(X; A; Q) = Q{Q,iX] A; Qo), and we may apply the results provided QiX (in- 
stead of X) satisfies the conditions. 

Theorem 2.1. (BG (1997a, Theorem 1.3)) Let 6 = 1/300, = 1^,3 = 13 and 
13 < c/ < oo. Let Co be an arbitrary positive absolute constant. Assume that condition 
Mqip, 6, So, CqG /a) holds. Then we have: 

A^^Uc{Ut + Al){l + \\a/af) (2.20) 

and 

A(;_Uc(n^ + A:)(l + ||a/af) (2.21) 
with C = cp"^ + c{cr/6s)^, where Of > 62 > . . . are the eigenvalues of (CQ)^. 

Unfortunately, we cannot apply Theorem 12.11 for d = 5, 6, . . . , 12. Moreover, the 
quantity C depends on p which is exponentially small with respect to eigenvalues of C. 

In Gotze and Zaitsev (2010), the following analogue of Theorem 12.11 is proved with 
bounds for constants which are not optimal. 

Theorem 2.2. Let 6 = 1/300, = 1^, s = 5 and 5 < d < 00. Let cq be an arbitrary 
positive absolute constant. Assume that condition AfQ{p, 6, So, cqG/ a) holds. Then 

A^N^ < C{a,'N-'/' E llXnf + a,' N'^ E ||X°f ) (l + \\a/af), (2.22) 

and 

A^nI < C {afE llXnf + a,-^iV-i E ||X°f ) (l + \\a/af), (2.23) 
with G = CdP^^. 

Theorem 12.21 extends to the case d > 5 Theorem 1.5 of BG (1997a) which contains the 
corresponding bounds for d > 9. Unfortunately, in both papers, the quantity G depends 
on p which is exponentially small with respect to Ug/o"^ (in BG (1997a)) and to a^/a'^ 
(in Gotze and Zaitsev (2010)). Under some additional conditions, G may be estimated 
from above by q exp(ccr^cr^^) and by q exp(ccr^cr^^) respectively. In Gotze and Zaitsev 
(2008) we proved Theorem 12.21 in the case a = and hence. Theorem 11.11 

The main result of the paper is Theorem 12.31 It is valid for 5 < c? < 00 in finite- 
dimensional spaces M.'^ only. However, the bounds of Theorem 12.31 depend on the smallest 
(jj's. This makes them unstable if one or more of coordinates of X degenerates. In our 
finite dimensional results. Theorems 12 . 3[ [2n and [2^8} we always assume that the covariance 
operator C is non-degenerate. 

Theorem 2.3. Let = I^^, 5 < d < 00. Then we have: 

<C(n° + A°)(l + ||a/(Tf ), (2.24) 

and 

Aj^Jn < C(n° + A°)(l + \\a/af), (2.25) 
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With C = Cd(T^ (detC)-i/2 _ 

It is easy to see that, according to (12.41) and (12.61) . 

n° + A°<E||C-i/2^f+V(rf^'+'^/'iV(^+'^/'), for 0<6<1, (2.26) 

and 

+ <E\\C-^/^Xf+^/{S^+^^/^N^/^), for 0<(5<2. (2.27) 
Therefore, Theorem 12.31 imphes the following Corollary 12.41 

Corollary 2.4. Let = 1^, 5 < < oo. Then we have: 

Ai;^ + lla/af )E||C-^/'Xf+V^^'^^^^', for < S < 1, (2.28) 

and 

Aj^Jn <dC(l + ||a/af )E||C-^/2Xf+VAr^/2, for 0<6<2, (2.29) 
with C = a'^ (det C)^^^"^ . In particular, 

max{AP,A^^]^} C ( 1 + ||a/af ) E ||C-^/2 ^|| V^V- (2.30) 

Theorem 12.11 and Corollary 12.41 yield Theorems 11.11 and II. 2[ using that -^0(2^) = 0, 
E xf < and + A| < n§ + < /3/(a^iV). 

Comparing Theorem 12.31 and Corollary 12 .41 with Theorem l2.3[ we see that the constants 
in Theorem 12.31 and Corollary 12 . 41 are written explicitly in terms of moment characteristics 
of C{X). In the case of non-positive definite quadratic forms Q such kind of estimates 
was unknown. 

If, in the conditions of Theorem 12. 3^ the distribution of X is symmetric or a = 0, then 
the Edgeworth corrections Ea{x) and El^{x) vanish and 

A^) = A^^]^ < C (n° + A°) (1 + \\a/af), C = caa" (det C)-'/^ (2.31) 

The corresponding inequality from Theorem 1.4 of BG (1997a) yields in the case s = 9 
and 9 < d < 00 under the condition Afqi^p, 6, So, coG/a) with 6 = 1/300 the bound 

Ai;) <C(n^ + A:) (1 + lla/af ), C = cp-\ (2.32) 

It is clear that sometimes the bound fl2.32p may be sharper than (I2.3ip . but unfortunately, 
it depends on p which is usually exponentially small with respect to crg/a"^. 

One can find more precise estimates of constants in the case of d-dimensional balls with 
d > 12 m the papers of Nagaev and Chebotarev (1999), (2005), Gotze and Ulyanov (2000), 
and Bogatyrev, Gotze and Ulyanov (2006). In this case Q = 1;^. See also Gotze and 
Ulyanov (2000) for lower bounds for A^^^ under different conditions on a and C{X). In the 
papers mentioned above, the authors have used the aproach of BG (1997a) and obtained 
bounds with constants depending on s < d largest eigenvalues erf > cr 2 >■■■> cr"^ of the 
covariance operator C (see Nagaev and Chebotarev (1999), (2005), with d > s = 13, and 
Gotze and Ulyanov (2000), and Bogatyrev, Gotze and Ulyanov (2006), with d > s = 12). 
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It should be mentioned, that, in a particular case, where Q = 1^ and d > 12, these results 
may be sharper than fl2.24p . for some covariance operators C. 

Thus, we see that the statement of Theorem 12.31 is especially interesting for d = 
5,6, ... , 11. It is new even in the case of d- dimensional balls. It is plausible that the 
bounds for constants in Theorem 12 . 31 could be also improved for balls with d > 5, especially 
in the case where d is large. It seems however that this is impossible in the case of 
general Q even if = 1^. For example, we can consider the operator Q such that 
Qcj = Cd-j+i, where Cej = crjCj, j = 1,2,. . . ,d, are eigenvectors of C. Following the 
proof of Theorem l2.3[ we see that the bounds for the modulus of the characteristic function 
I "ifbit) I = I E e{ t Q [Zj^ ~ ^] } I behave as the bounds for the modulus of the characteristic 
function | E e\^tId[Z]\f — b] } | but with eigenvalues of the covariance operator aiad, o'2crd^i, 
(T3(Td_2, . . . which can be essentially smaller than o"^ > (t| > cTg > ■ ■ ■ . Therefore, it is 
natural that the bounds for constants in Theorem 12.31 depends on the smallest eigenvalues 
of the covariance operator C 

Note that, in the proof of Theorem [O in BG (1997a), inequalities fl2:20D and fl2:2T]) 
were derived for the Edgeworth correction Ea{x) defined by (I2.13p . However, from The- 
orems 12.11 and 12.21 or 12.31 it follows that, at least for 13 < d < oo, definitions fl2.13p 
and (I2.14P determine the same function Ea{x). Indeed, both functions may be repre- 
sented as N~^^'^ K{x), where K{x) are some functions of bounded variation which are 
independent of A^. Furthermore, inequalities fl2.20p and fl2.24p provide both bounds of 
order 0{N~^). This is possible if the Edgeworth corrections Ea{x) are the same in these 
inequalities. 

On the other hand, it is proved (for d >9) that definition 02.130 determine a function 
of bounded variation (see BG (1997a, Lemma 5.7)), while definition (12.140 has no sense 
for d = oo. 

Introduce the concentration function 

Q{X; \) = Q{X; X;Q) = sup P{x < Q[X - a] < x + \}, for A > 0. (2.33) 

a, xeM<* 

It should be mentioned that the supremum in (I2.33P is taken not only over all x, but 
over all x and a G M'^. Usually, one defines the concentration function of the random 
variable Q[X — a] taking the supremum over all x G M*^ only. Note that, evidently, 
Q{X + Y; X) < Q{X; A), for any Y which is independent of X. 

The following Theorems 12.51 and Theorem 12.61 are Theorems 1.5 and 2.1 from Gotze 
and Zaitsev (2010). 

Theorem 2.5. Let Q"^ = Id, 5 < s < d < oo, s < oo and < 5 < l/(5s). Then we have: 
{i) If condition A/q(p, 6, So, X) is fulfilled with some p > 0, then 

Q{Zn; A) (pN)-^ max{l; A}, for A > 0. (2.34) 

(ii) If, for some m, condition J\fQ(p,S, So, m^^^'^ Z^) is fulfilled, then 

Q{Zn; A) ipN)~^ max{m; A}, for A > 0. (2.35) 
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Theorem 2.6. Let = 1^ and 5 < d < oo. Let cq be an arbitrary positive absolute 
constant. Assume condition //qIp, 6, So, cqG /a) to be satisfied with s = 5 and 6 = 1/200. 
Then 

Q{Zn; X) '^p''^max{Ul + Al; Xa-^N-^}, forX>0. (2.36) 
In particular, Q{Zn] A) <^p^'^N^^ max{ /Ja^"^; Acr^^}. 

Theorems 12.51 and Theorem 12.61 extend to the case 5 < < oo Theorems 1.6 and 2.1 
of BG (1997a) which were proved for 9 < d < oo. 

We say that a random vector Y is concentrated in L C M'' if Y*{Y G L} = 1. In 
BG (1997a, item {Hi) of Theorem 1.6) it was shown that if X is not concentrated in a 
proper closed linear subspace of M'^, 1 < d < oo, then, for any 5 > and S there exists a 
natural number m such that the condition Mq{p, 5, S, m~^^'^ Zm) holds with some p > 0. 

In this paper, we shall prove the following Theorems 12.71 and Theorem 12.81 

Theorem 2.7. Let Q'^ = I^, 5 < s = d < oo and < 5 < l/(5s). Then we have: 
(i) If condition A/'(p, 5, So, C"^/^ X) is fulfilled with some p > 0, then 

Q{Zn; A) <d (pN)-^ max{ 1; A } a'^ (det C)~l/^ for A > 0. (2.37) 
(m) If , for some m, condition J\f(p,S, So, m^^^'^C^^^'^ Z^) is fulfilled, then 

Q{Zn; A) <d (pN)-^ max{m; A (j-^ } (x'^ (det C)-^/^ for A > 0. (2.38) 
Theorem 2.8. Assume that 5 < d < oo and that = 1^. Then 

Q(Z^; A) <dmax{n° + A°; A(T"^X"^}(T'^(detC)"^/^ forX>0. (2.39) 
In particular, Q{Zn] A) <d N'^ max{ E ||C-i/2 j^f ; Xa-'^] {deiC)-^/^ . 

Theorem 12.71 and Theorem 12.81 yield more explicit versions of Theorems 12.51 and The- 
orem 12.61 as well as Theorem 12.31 is in a sense a more explicit version of Theorem 12. 2[ 
We should mention that Theorems 12.21 12.51 and 12.61 do not follow from Theorems 12. 3^ 12.71 
and 12.81 For the proofs of these theorems we refer the reader to the paper of Gotze and 
Zaitsev (2010) and to the preprint of Gotze and Zaitsev (2009) which are available in 
internet. 

For example, the bounds in Theorems 12. 2[ [275] and [2y6] may be sharper than those from 
Theorems 12.31 12.71 and 12.81 in a particular case, where Q = 1^ and x^^ a. Under some 
additional conditions, a'^ (det C)~^/^ is replaced by exp(ccr~^crf ) 1. On the other hand, 
0"'' (detC)^^/^ provides a power- type dependence on eigenvalues of C and the results are 
valid for Q which might be not positive definite. 

In Theorems 12.31 and Theorem l2.8l we do not assume the fulfilmemt of conditions M{-) 
or A/'q( ■ ). In the proofs, we shall use, however, that, for an arbitrary absolute positive con- 
stant Co and any positive quantity q depending on d only, condition ^{p, 5, So, Cq C^^/^ G) 
is fulfilled with s = d, 5 = Cd and p x^ 1, for any orthonormal system So- 

Similarly to BG (1997a), in Section [3], we prove bounds for concentration functions. 
The proof is technically simpler as that of Theorem 12. 3[ but it shows how to apply the 
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principal ideas. This proof repeats almost literally the corresponding proof of BG (1997a). 
The only difference consists in the use of new Lemma 18.31 which allows us to estimate 
characteristic functions of quadratic forms for relatively large values of argument t. In 
Sections m and Theorem 12.31 is proved. We shall replace Lemma 9.4 of BG (1997a) by 
its improvement, Lemmas 15.11 Another difference is in another choice of k in fl5.3ip and 
(I5.32P in comparison with that in BG (1997a). In Sections |6HH] we prove estimates for 
characteristic functions which were discussed in Section [TJ 

Acknowledgment We would like to thank V.V. Ulyanov for helpful discussions. 



3. Proofs of bounds for concentration functions 

Proof of TheoremsWTA and \2.8l Below we shall prove the assertions fl2.37p : fl2.37p =^ 
fl2.38p and fl2.38p =^ fl2.39p . The proof repeats almost literally the corresponding proof 
of BG (1997a). It is given here for the sake of completeness. The only essential difference 
is in the use of Lemma \K3\ in the proof of Lemma [3711 We have also to replace everywhere 
9 by 5 and o by □. □ 

For <tQ <T and b eM.'^, define the integrals 



\t\ 



Io = I \ ^b{t) \dt, h= I $b(t) 

'to<\t\<T 



where 

$fe(t) =Ee{tQ[Z;v-&]} (3.1) 
denotes the Fourier-Stieltjes transform of the distribution function \E'b of Q\Z]^ — h] (see 
fim]) and fITTTp i Note that \ ^b{-t) \ = \ ^b{t)\- 

Lemma 3.1. Assume condition J\f(p, S, So, C^^^^ X) with some < 6 < l/(5s) and 
5 < s = d < oo. Let cr^ = 1 and 

to = c,{s)a^\pN)-'+'/', C2{s)a^^ <T< c^is) a^' (3.2) 

with some positive constants Cj{s), 1 < j < 3. Then 

lo (det C)-^/2 ipN)-\ h <s (det C)"^/' ipN)-\ (3.3) 

Proof. Note that the condition cr^ = 1 implies that 

T X, al X, (T^ = 1 and det C 1. (3.4) 

Denote k = pN. Without loss of generality we assume that k > Cs, for a sufficiently large 
quantity Cs depending on s only. Indeed, if k < Cg, then one can prove (13. 3p using (13. 4p 
and |\E'f,| < 1. Choosing Cs to be large enough, we ensure that k > Cg implies 1/k < to < T. 
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Lemma 18.31 and (13. 4p imply now that 

/ h^^^W — «. ^ ^ , (3.5) 

for any 04(5) depending on s only. Inequalities (13. 4p and (13. 5p imply (13. 3p for Ji. 

Let us prove inequality (13. 3p for Jq. By (13. 4p and by Lemma WT\ for any 7 > and 
any fixed t G M satisfying k^^"^ \t\ < 05(3), where 05(3) is an arbitrary quantity depending 
on s only, we have (taking into account that |\E'fe| < 1) 

I ^b{t) I min{ 1; + k-''^ \t\-''^ (det C)~^'^], k = pN. (3.6) 

Furthermore, choosing an appropriate 7 and using (13. 40 - 03. 6p . we obtain 

pl/k -, poo J, 1 

(det Cf'^ Jo «. I dt+-+ j^^^ ^ «, - , (3.7) 

proving (13.30 for /q. □ 

Proof of (I2.37p . Let o"^ = 1. Using a well-known inequality for concentration functions 
(see, for example, Petrov (1975, Lemma 3 of Ch. 3)), we have 



:{A;1} 

Jo 



<5(^jv; A) < 4 sup maxj A; 1 1 / \'^b{t)\dt. (3.8) 



To estimate the integral in (13. 8 p we shall apply Lemma 13.11 which implies that 

Q{Zn; A) <rfmax{A; l}{pN)-\detCy^/^, (3.9) 

proving (12.370 in the case cr^ = 1. If cr^ 7^ 1, we obtain (I2.37P applying (13. 9 p to ZN/a. □ 

Proof of (12.371) =^ (12.381) . Without loss of generality we can assume that N/m > 2. 
Let Yi, • • • be independent copies of m-i/2^„. Denote Wk = Yi + ■ ■ ■ + Yk. Then 
C{Zn) = C{y/mWk + y) with k = \N/m\ and with some y independent of Wk- Therefore, 
Q{Zn] A) < QiWk] X/m). In order to estimate QiWk, A/m) we apply (12.370 replacing Zn 
by Wk- We have 

Q{Wk;X/m) (pA;)-^max{l; Aff-Vm} a'^ (detC)"^/^ 

{pN)-^m8ix{m; Xa-^ja"^ {det Cy^/\ □ (3.10) 



Recall that truncated random vectors and their moments are defined by (I2.3p - (l2.6p 
and that C = cov X = cov G. 

Lemma 3.2. The random vectors X^, X^ satisfy 

{Cx,x) = (covX°x,x) + E(Xn,x)^ + (EX°,x)^ 
There exist independent centered Gaussian vectors and W such that 

C{G) = £(G, + W) 
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and 

2covG', = 2covX° = covXS, {coyWx.x) = E (Xn, a;)^ + (EX°, x)^. 
Furthermore, 

E ||C-^/2 ^||2 = ^ = E \\C-^'^ + E ||C-^/2 p^||2 

an(i E||C-i/2p^f < 2rfn°. 

We omit the simple proof of this lemma (see BG (1997a, Lemma 2.4) for the same 
statement with o instead of □). Lemma [3.21 allows us to define the vector X' by (12.101) . 

Recall that and Z]^ denote sums of independent copies of and X^ respec- 
tively. 

Lemma 3.3. Let e > 0. There exist absolute positive constants c and C\ such that the 
condition < Cip6'^/{de^) implies that 

N{p, 6, S, e G) =^ Ar(p/4, 45, 5, e {2m)-^/^ 
form> ce^d^N A^/{p6'^). 

Lemmas 13. 21 and 13.31 are in fact the statements of Lemmas 2.4 and 2.5 from BG (1997a) 
applied to the vectors C^^^^X instead of the vectors X. We use in this connection 
equalities ([22D, (I23D and (^M> replacing in the formulation A^, B^, G, Z^, . . . by d, 
A°, B°, G, Z^, ... respectively. 

Proof of (12.381) =^ (I2.39p . By a standard truncation argument, we have 

\P{Zn e A] -P{^^ e A} I < XP{ \\C-^/^X\\ > VdN} < B°, (3.11) 
for any Borel set A, and 

g(Z^,A)<B° + Q(Z°,A). (3.12) 

Recall that we are proving (I2.39P assuming that 5 < < oo. It is easy to see that, for an 
arbitrary absolute positive constant Cq, condition N{p,S,So,CoC^^^'^ G) with 

s = d, (5= 1/(20 s), p-dl (3.13) 

is in fact fulfilled automatically for any orthonormal system So, since the vector C^^^^ G 
has standard Gaussian distribution in M.'^ and P| ||coC^^/^G — e|| < 5} = c{d) for any 

vector e eW^ with ||e|| = 1. Clearly, 45 = l/(5s). Write K = e/V2 with e = cq. Then, 
by f l3.13p and Lemma 13. 3[ we have 

Af{p, 6, So, ^) ^ ^-1/2 ^c-i/2 (3^14) 

provided that 

n^<ci{d), m>C2{d)NA^. (3.15) 

Without loss of generality we may assume that B° < ci{d), since otherwise the result 
follows easily from the trivial inequahty Q{Zn', A) < 1. 
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The non-degeneracy condition fl3.14p for KZ'^^ allows to apply fl2.38p of Theorem 12.71 
and, using f l3.13p . we obtain 

g(Z^,A) = Q{KZ^,K^X) <rfAr-imax{m; fsT^ A/KV^} (det C)"^/^ (3.16) 

for any m such that f l3.15p is fulfilled. Choosing the minimal m in f l3.15p . we obtain 

g(Z°, A) <rf max{A°; X/ia^N)} a'' (det C)"^/^ (3.17) 

Combining the estimates fl3.12p and fl3.17p . we conclude the proof. □ 



4. Auxiliary lemmas 

In Sections H] and E] we shall prove Theorem 12.31 Therefore, we shall assume that 
its conditions are satisfied. We consider the case d < oo assuming that the following 
conditions are satisfied: 

q2 = 1^^ = 1, d>5, b = ViVa. (4.1) 

Moreover, it is easy to see that, for any absolute positive constant cq and for any or- 
thonormal system So = {ei, . . . , e^} C M'^ , condition 

N{p,S,So,coC-^^'^G) with j9 1, 5 < s = d < oo, 5 = 1/(20 s) (4.2) 

is in fact fulfilled automatically since the vector C^^/^ G has standard Gaussian distribu- 
tion in and, therefore, P{ ||co C'^/^ ^ - e|| < 5} = P{ || C'^/^ ^ - Cq ^ e|| < c^^ 6} = 
c{d) for any vector e G M"' with ||e|| = 1. 

Notice that the assumption o"^ = 1 does not restrict generality since from Theorem 12.31 
with 0"^ = 1 we can derive the general result replacing X, G by X/a, G/a, etc. Other 
assumptions in (14.10 are included as conditions in Theorem 12.31 Section H] is devoted to 
some auxiliary lemmas which are similar to corresponding lemmas of BG (1997a). 

In several places, the proof of Theorem 12.31 repeats almost literally the proof of The- 
orem 1.5 in BG (1997a). Note, however, that we shall use truncated vectors X?, while 
in BG (1997a) the vectors Xj were involved. We start with an application of the Fourier 

transform to the functions ^i, and $5, where b = a//V a. We shall estimate integrals 
over the Fourier transforms using results of Sections [31 EHHl and some technical lemmas of 
BG (1997a). We shall also apply some methods of estimation of the rate of approximation 
in the CLT in multidimensional spaces (cf., e.g., Bhattacharya and Rao (1986)). 

Below we shall use the following formula for the Fourier inversion (see, for example, 
BG (1997a)). A smoothing inequality of Prawitz (1972) implies (see BG (1996, Section 4)) 
that 

Fix) = 1 + ^ V. P. 1^^^^ e{ -xt}F{t) ^+R, (4.3) 
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for any K > and any distribution function F with characteristic function F (see fl2.11l) ). 
where 

\R\<-^[ \m\dt. (4.4) 

'\t\<K 



Here V. P. / f{t) dt = Utcis^o i|t|>e denotes the Principal Value of the integral. 

Recall that the random vectors X' are defined in (12. 4p and (12.101) and Z^, Z'^ are 
sums of their independent copies. Note that the Gaussian vector W involved in (I2.10p 
is independent of all other vectors and have properties described in Lemma 13.21 Write 
and \E''{, for the distribution function of Q[2'^ — h] and Q[^^ — h] respectively. For 
< k < N introduce the distribution function 

^f\x) = P{ Q [Gi + ■ ■ ■ + + X^+i + ■ ■ ■ + - 6] < X } . (4.5) 

Notice that = %, = $6. 

The proof of the following lemma repeats the proof of Lemma 3.1 of BG (1997a). The 
difference is that here we use the truncated vectors instead of XJ. 

Lemma 4.1. Let Cd be a quantity depending on d only. There exist positive quantities ci{d) 
and C2{d) depending on d only such that the following statement is valid. Let < Ci{d)p 
and let an integer 1 < m < N satisfy m > C2{d) NA^/p, Write 

K = ell {2m), t, = ca{pN/m)~^+^'''. 

Let F denote any of the functions \E'^, \E'^, '^^^^ or Then we have 

F{x) = 1 + ^ V. P. /" e{-xtK} F{tK) ^ + Ri, (4.6) 

with |i?i| <d (piV)-im (detC) -1/2 _ 

Proof. We shall assume that {pN)~^m < c^ld) with sufficiently small c^ld) since 
otherwise the statement of Lemma 14.11 is trivial (see (13. 4p . (14. 3 p and (I4.4p ). Let us 
prove (14. 6p . We shall combine (14. 3 p and Lemma 13.11 Changing the variable t = tK 
in formula (14. 3p . we obtain 

F{x) = 1 + ^ V. P. /" e{-xtK} F{tK) ^ + R, (4.7) 

>^|t|<i 

where 

\R\< [ \F{tK)\dt. (4.8) 

J\t\<i 

Notice that "^[^ and are distribution functions of random variables which 
may be written in the form: 

+ T] , v = G^ + --- + Gk + X^_^, + --- + X^, 
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with some k, < k < N , and some random vector T which is independent of Xj and Gj, 
for all j. Let us consider separately two possible cases: k > N/2 and k < N/2. 

The case k < N/2. Let Y denote a sum of m independent copies ofK^/^X^. Let 
Yi, Y2, ... be independent copies of Y. Then we have 

dK^^W) = C(Yi + --- + Yi + Ti) (4.9) 

with / = \N/{2m)'] and some random Ti independent of Yi,...,Yj. By ( 14.21) and by 
Lemma 13.31 we have 

Afip, 6, S, coC~^/2 ^ Afip/A, 46, S, C^^/^ Y) (4.10) 

provided that 

n°<p/c/3 and m:^ d^NA^/p. (4.11) 

The inequalities in ( 14. lip follow from conditions of Lemma 14.11 if we choose some suffi- 
ciently small (resp. large) ci{d) (resp. C2{d)). Due to ( lO) . f Ol) . ( 149|) and f l4.10p . we 
can apply Lemma [3.11 in order to estimate the integrals in (14. 7p and (14. 8p . Replacing in 
Lemma [3. II X by Y and by /, we obtain (14. 6 p in the case k < N/2. 

The case k > N/2. We can argue as in the previous case defining now y as a sum of 
m independent copies of K^^'^G. Condition A/'(p/4, 45, 5o, C~^/^ F) is satisfied by (14.21) . 
since now C{Y) = C{cqG). □ 

Following BG (1997a), introduce the upper bound A^,£(X),£(G)) for the char- 
acteristic function of quadratic forms (cf. Bentkus (1984) and Bentkus, Gotze and Zitikis 
(1993)). We define A^, £(X), £(G)) = x(t; A^, £(X)) + x(t; A^, £(G)), where 

x(t;X,£(X)) = sup |E e{tQ[Zj-] + |, Zj = Xi + ■ ■ ■ + Xj, (4.12) 

with j = ^ (N — 2)/14] . In the sequel, we shall use that 

x(t; N, C{X'), C{G)) < Kit- N, £(X°), C{G)). (4.13) 
For the proof, it suffices to note that X' = X^ — E X^ + W and W is independent of X^. 

Lemma 4.2. Let the conditions of Lemma \4.1\ be satisfied. Then 

[ {\t\KyK{tK-,N,C{X^),C{G)) ^ 

J\t\<ti I*-! 

{{Npy, for 0<a< d/2, 

{Npy {l + \\og{Np/m)\), for a = d/2, (4.14) 
{Npy (1 + (Xp/m)(2"-^)/'^), for a > d/2. 

Lemma IT2] is a generalization of Lemma 3.2 from BG (1997a) which contains the same 
bound for < a < d/2. In this paper, we have to estimate the left hand side of (14.141) in 

the case d/2 < a too. 
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Proof. We shall assume again that {pN)^^m < €^{(1) with sufficiently small c^ld) since 
otherwise fl4.14p is an easy consequence of |x | < 1. 

Note that \ B e{tQ[Zj] + {x, Zj)}\ = \ E e{tQ[Zj - y]}\ with y = -Qx/2. By dMD 

and flCTD . the condition Af{p/A,46,So, K^^^C-^^^ Z^^) is fulffiled. Therefore, collecting 
independent copies of K^/^X° in groups as in ( 14. 9p . we can apply Lemma [8.11 By f l3.4p . 
(14 .2^ and this lemma, for any 7 > and \t\ < ti, 

>i{tK;N,C{X°)) -^^4 {pN/ my ^ + mm{l; {Np/m)-'^'^ \t\-'^'^ (detC)-^/2j^ 

We have used that = 1 implies cr^ 1. A similar upper bound is valid for the quantity 
>f{tK]N,C{G)) (cf. the proof of gl]) for k > N/2). Thus, we get, for any 7 > and 

\t\ < h, 

x{tK-N,C{X^),C{G)) (j9iV/m)-^ + min{l; (detC)-^/^ (^ni/ {\t\pN)Y'''] . 
Integrating this bound (cf. the estimation of Ji in Lemma [3. II) . we obtain (I4.14p . □ 



5. Proof of Theorem 12.31 

To simplify notation, in Section [5] we write 11 = 11^ and A = A4 . The assumption 
(T^ = 1 and equalities ^WC'^'"^ XW^ = ([MD and ([2SD imply 

n + AA^>l, n + A<l, (t| < 1, detC<l. (5.1) 

Recall that aJ^^ and functions ^b, $6 and 9^ are defined by ^HM and fl2Tn) -f l239l ). 
Note now that e^(x) = E^{x/N) and, according to (KW\ . 

A^^^ < A^;)^ + sup \Q,{x) - e°(a:) |, (5.2) 

where b = a/ZV a and 

A^^^n = sup I ^fe(x) - $5(3;) - Q°ix)\. (5.3) 

xGR 

Let us verify that 

sup |eb(a;)-e°(x)| «,n°. (5.4) 

x£R 

To this end we shall apply representation f l2.14p -( l2.15p of the Edgeworth correction as a 
signed measure and estimate the variation of that measure. Indeed, using ( I2.14p - (l2.15p . 
we have 

sup |e6(x)-0^(a;)| < AT-i/'j, I = [ \Ep'"{x)X^ -Ep"'{x)X''^\dx. (5.5) 

xm jRd 

By the explicit formula (12.161) . the function u 1— )■ p"'{x)u^ is a 3-linear form in the vari- 
able u. Therefore, using X = X^ + Xo and ||Xn|| = 0, we have p"'{x)X^ — 
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p"'(a;)X°' = p"'{x)X^, and 

iV-i/2/< 3^=^/2 nn f [\\C-'/^x\\ + \\C-'/^xf)p{x)dx = CdU'^. (5.6) 

Inequalities (15 .Sp and (15. 6p imply now (15. 4p . 

To prove the statement of Theorem 12.31 we have to derive that 

A% «, (n + A)(l + \\a\\Y (det C)-'/'. (5.7) 

While proving (15.71) we assume that 

n < Cd, and A < q, (5.8) 

with a sufficiently small positive constant Cd depending on d only. These assumptions do 
not restrict generality. Indeed, we have | \I^b(x) — $b(a;) | < 1. If conditions (15. 8p do not 
hold, then the estimate 

sup |e°(a;)| «,Ar-i/2E||C-i/2x°f «,Ai/2 (5.9) 

immediately implies (15. 7p . In order to prove (15. 9p we can use (12. 6 p and representation 
(I2.14p - (l2.15p of the Edgeworth correction. Estimating the variation of that measure and 
using 

E ||C-^/2 < E \\C-^/^Xf = d, (5.10) 

(E X°f )2 < E X°f E , (5.11) 

we obtain (15. 9p . 
It is clear that 

A^Ja < sup ( I v^,(x) - %ix) I + I e°(a;) - G^x) | + | %ix) - <^,{x) - Q'.ix) \ ) . (5.12) 
Similarly to (15.51) . we have 
sup |e^(x) - e',{x) I < J, J = [ \Ep"'{x)X^^ - Ep'"{x)X'^\dx. (5.13) 

Recall that vector X' is defined in (IXTOj) . By Lemma [321 we have E HC^^/^ VTf < 2c/n 
(hence, E ||C~^/^ 11^/^, for < g < 2). Moreover, representing W as a sum of 

a large number of i.i.d. Gaussian summands and using the Rosenthal inequality (see BG 
(1997a, inequality (1.24)), we conclude that 

EllC-^/'iyf (EllC-^/'PFll^)"/' <,,^n'?/2, g>0. (5.14) 

Furthermore, according to ( 12. 4p . ( 12. 6p and (15. Sp . 

Hence, by ([M]), ([210]), ([EI]), ([EH and ([EIS]), 

E llX'f < ^ = E ^/||4 ^ (5.16) 
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Using dsn, dSID, O, dSini) and (^-(^, we get 

N-'/^J ni/2(iV-i/2n + AV2) /" (||C-i/2a;|| + ||C-i/2xf)p(x)tia: 

J«.<i 

<d n + A. (5.17) 
Thus, according to fl5.13p and fl5.17p . 

sup |e^(x) -e;(x)| <dn + A. (5.18) 



The same approach is apphcable for the estimation of | ©j, | • Using f l2.10p . ( I2.14p -( l2.16p . 
CT. dSmi), (imi), dsn and (15151), we get 

sup I e;(x) I < f I E/'(x)X'^ I rfx 

«, Ai/2 + Ar-V2n3/2. (5.19) 

Let us prove that 

sup I'^bix) -%{x) \ < (detC)~i/2p-2(n + A)(l+ ||af). (5.20) 

Using truncation (see (13. lip ), we have — < 11, and 

sup I ^b(x) - %{x) I < n + sup I *^(x) - ^fe(a;) I- (5.21) 

In order to estimate |\&^ — ^fel, we shall apply Lemmas 14.11 and 14.21 The number m in 
these Lemmas exists and NA/p 1, as it follows from (15. ip and (15. 8p . Let us choose the 
minimal m, that is, m NA/p. Then {pN)~^m A/p^ and m/A^ A/p. Therefore, 
using Lemma [4.11 we have 

sup I vl/°(x) - %{x) I «, A (det C)-i/2 + / I $°(r) - $;(r) | ^ , r = 

(5.22) 

We shall prove that 

|$°(r) - %{t) I <rf xn |r| Ar(l + \r\N){l + ||a|p), (5.23) 

with X = x(r; A^, £(X'-')). Combining (I5.2ip -( l5.23p . using r = ti^ and integrating in- 
equality (I5.23P with the help of Lemma 14.21 we derive (I5.20p . 

Let us prove (15:231) . Recall that X' = X° - EX° + PF, where W denotes a centered 
Gaussian random vector which is independent of all other random vectors and such that 
covX' = covG (see Lemma [3.20 . Writing D = — EZ^ — b, we have 

Z°-6 = D + EZ°, C{Z'j^-b) = C{D + VnW), 

and 

|$°(r)-$Ur)| < |/i(r)| + |/2(r)| (5.24) 
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with 

/i(r) = E e{TQ[D + VnW] } - E e{rQ[D] }, 

/2(t) = E e{ r Q + E } - E e{ r Q [D] } . 

Now we have to prove that both | /i(t) | and | /2(t) | may be estimated by the right hand 
side of (IE23]). 

Let us consider /i. We can write Q[D + y/NW] = Q[D] + A + B with A = 
2\/N (Q-D, W) and B = NQ[W]. Taylor's expansions of the exponent in fl5.25p in pow- 
ers of irB and irA with remainders 0{tB) and 0{t'^ A^) respectively imply (recall that 
E = and = 1^) 

|/i(r)| < x|r|ArE||iyf + xr2ArE||Vrf E||L)||2, (5.26) 

where x = x(r; A^, £(X°)). The estimation of the remainders of these expansions is 
based on the splitting and conditioning techniques described in Section 9 of BG (1997a), 
see also Bentkus, Gotze and Zaitsev (1997). Using the relations cr^ = 1, E||IVp ^ 
E ||C-i/2y^f <rf n and E \\Df < A^(l + \\af), we derive from (E^S]) that 

|/i(r)| <,xn|r|iV(l + |r|iV))(l + ||af). (5.27) 

Note that EZ^ = A^EX° = -N^Xu- Expanding the exponent e{rQ[L) + EZ^]}, 
using ( I5.15P and proceeding similarly to the proof of fl5.27p . we obtain 

|/2(r)| xn|r|iV(l + ||a||). (5.28) 

Inequahties (EMD, and fl^:^ imply now ([523]) • 

It remains to estimate | \E'^ — — 6J, |. Recall that the distribution functions '^f^{x)^ 
for < / < A^, are defined in 

Fix an integer k, 1 < k < N. Clearly, we have 

sup I %ix) - - e;(x) I < Ji + /2 + /3, (5.29) 

where 

/i = sup l^f^(x) - <i>bix) -{N-k) Q'f,{x)/N\, (5.30) 

/2 = sup|^;(x)-^f)(x)|, (5.31) 

and 

h = sup kN-^ \Q'b{x)\. (5.32) 

xGR 

Let estimate Ji. Define the distributions 

N 

fi{A) = P{Uk + 5^ e Vna}, fio{A) = p{UNe Vna} = P{G e A}, 
j=k+i 

(5.33) 
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where Ui = Gi + ■ — \- Gi. Introduce the measure x' replacing X by X' in f l2.15p . For the 
Borel sets A C M'^ define the Edgeworth correction (to the distribution /i) as 

ljf\A) = {N- k)N-^/WA)/Q. (5.34) 
Introduce the signed measure 

p = jj, — jiQ — IJif\ (5.35) 
It is easy to see that a re- normalization of random vectors implies (see relations f l2.14p . 
^M- drii, ill) and (EJl-dESS])) 

\^f\x) -%{x) - {N -k)Q'^^{x)/N\ = iy{{u e R"^ : Q[u - a] < x/N}) 

< 6n = sup |z/(A)|. (5.36) 



Lemma 5.1. Assume that d < oo and 1 < k < N . Then there exists a c{d) depending 
on d only and such that 5n defined in fl5.36p satisfies the inequality 

5N<^d + exp{ -c(rf) k/(3 } (5.37) 

mt/i ^ = E||C-i/2X'f . 

An outline of the proof. We repeat and slightly improve the proof of Lemma 9.4 in 
BG (1997a) (cf. the proof of Lemma 2.5 in EG (1996)). We shall prove (I5.37P assuming 
that covX = covX' = covG = 1^, Applying it to C'^/"^ X' and C'^/'^G, we obtain (KST} 
in general case. 

While proving (15.371) we assume that P/N < Cd and N >l/cd with a sufficiently small 
positive constant q. Otherwise fl5.37p follows from the obvious bounds j3 > a'^ = d"^ and 

6n <^dl + {P/Ny/^ I \\xfp{x)dx<^dl + (P/NY/\ 

Set n = N — k. Denoting by and U'j sums of j independent copies of X' and G' 
respectively, introduce the multidimensional characteristic functions 

g{t) = 'Ee{{N-^'H,G)], /i(t) = E e{ (AT-^/^t, X') }, (5.38) 

f{t) = E e{ {N-'IH, Z;) } = h-{t), fo{t) = E e{ {N-'^H, K) } = g^it), (5.39) 

/i(t) = nm{t) /o(t), where m{t) = E {tt,X')\ (5.40) 

Ht) = (fit) - fo{t) - /i(t)) g{pt), p' = k. (5.41) 

It is easy to see that 

dit) = / e{{t,x)}u{dx). (5.42) 
Using the truncation, we obtain 

E||z;/yiV|r <^,rf 1, 7>0, 1</<X. (5.43) 
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By an extension of the proof of Lemma 11.6 in Bhattacharya and Rao (1986), see also 
the proof of Lemma 2.5 in BG (1996), we obtain 

6n <d max / I d^u^t) I dt. (5.44) 

\a\<2d Jj^d 

Here |a| = |ai| + ■ ■ ■ + la^l, a = (ai, . . . , ad), G Z, aj > 0. In order to derive (15.37^ 
from f l5.44p . it suffices to prove that, for \a\ < 2d, 

Id'^Ht)] gicipt), (5.45) 

|5"z^(t)| <d PN-\l + \\t\f) exp{-C2\\tf}, for \\tf <C3{d)N/ p. (5.46) 



Indeed, using f l5.45p and denoting T = \J c^{d)N/(5, we obtain 

/ 1 9"z>(t) I <d / g{cipt)dt<^d exp\^- j exp{-cl\\t\\y8}dt, 

\\t\\>T \\t\\>T K'^ 

(5.47) 

and it is easy to see that the right hand side of (15.471) is bounded from above by the 
second summand in the right hand side of fl5.37p . Similarly, using fl5.46p , we can integrate 
|9"P(t) I over < T, and the integral is bounded from above by Cdl3/N. 

In the proof of (15. 450 - 05. 47p we applied standard methods of estimation which are 
provided in Bhattacharya and Rao (1986). In particular, we used a Bergstrom type 
identity 

n— 1 n— 1 j — 1 

/ - /o - /i = Y.^h -g-m) W g^-^-^ + ^ m ^^(/i - g) g--'-\ (5.48) 

i=0 i=0 1=0 

relations flSlSHjl - flCT]) . 1 < A; < iV, | 9" exp{-C4 \\tf} \ exp{-C5 \\tf}, \fN jf^ >d 1 
and y'^'^ exp{— 1, for ?/ > 0. □ 

Applying (I5.30p . (I5.36P and Lemma [STTl we get 

^1 ^ + exp{-c(d)A;//3}. (5.49) 

For the estimation of we shall use Lemma 15.21 which is an easy consequence of BG 
(1997a, Lemma 9.3), fl4A3l) and flKITHl) . 

Lemma 5.2. We have 



|$;(t)-$['^(t)| <xt2/(/3 + |t|A^/3 + |t|A^yA^/3)(l + ||af), /or < / < A^, 
where Kit; N,L{X^),L{G)) {of. Km ). 
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As in the proof of fl5.22l) . applying Lemma l4.ll (choosing m A^(A + n)/p) and 
using ( I4.2p . we obtain 

/2«rf(A + n)(detC)-i/^+ [ \%{T)-¥,'\T)\dt/\t\, T = tK. 

J\t\<ti 

The existence of such an m is ensured by f l4.2p , (15.1 p and (15. 8 p , Applying Lemma 15.21 and 
replacing in that Lemma t by r, we have 

|$;(r) - $f^(r)| < >cT'^k{l3 + \T\N^+\T\N\[Nj){l+\\af). (5.50) 
Integrating with the help of Lemma [4.21 and using (14.20 . we obtain 



h <d (detC)-i/2(n + A + A;Ar-2(/3 + ^Ar/3)(i + (n + A)-i/'^)(l + ||af)). (5.51) 

Let us choose k N^/'^jf^'^ . Such k < N exists by ^ = 1, by (I5.16P and by 

assumption (15. 8p . Then (15.490 and (15.510 turn into 

and 

h «, (det C)-^/^ (n + A + [{^Y" + [-^f") (1 + (n + A)-V<^) (1 + llaf )). (5.53) 

Using (112D, (lElD, fl5:T6|) and fl533|) . we get 

/2 «d (det C)-i/2 (n + A + ^ (1 + ||af )). (5.54) 

Finally, by (EH]), (EIED, f lCTD and (^M>, 

h «d ^ (A^/2 + iV-V2n3/2) « A + n. (5.55) 

Inequalities ([EHD, dSig), flSTTG]) . flSTTSj) . fl5:20|) . flOO]) . fl532|) . f lCTD and fl535|) imply 
now (15. 7p (and, hence, (I2.25P ) by an application of 11 + A < 1. Note that, by (12. 6p . we 
have n < n!^. Together with (ED and (El), inequality (EZD yields ^I?2M . The statement 
of Theorem 12.31 is proved. □ 



6. From Probability to Number Theory 

In Section [6] we shall reduce the estimation of the integrals of the modulus of charac- 
teristic functions '^b{t) to the estimation the integrals of some theta-series. We shall use 
the following lemmas. 
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Lemma 6.1. (BG (1997a, Lemma 5.1)) Let L, C g M'^ and let Q : R'^ ^ M.'^ be a 
symmetric linear operator. Let Z, U, V and W denote independent random vectors taking 
values in M'^. Denote by 

P{x) = {Qx,x) + {L,x) + C, ,tgM^ 

a real-valued polynomial of second order. Then 

E e{tP{Z + U + V + W)} ^ < E e{2t{QZ,U)} + E e{2t{QZ,V)}. 

Let ei,e2,... denote i.i.d. symmetric Rademacher random variables. Let 6 > 0, 
S = {ei, . . . , e^} C M"^ and let D : M'' — )■ M*^ be a linear operator. Usually, we shall take 
D = C^"*^/^. We shall write C{Y) G r(5; ©, S) if a discrete random vector Y is distributed 
as £1^1 + ■ ■ ■ + EgZs , with some (non-random) Zj G M'^ such that \\B)Zj — ej\\ < 6, for all 
1 ^ J ^ -5. Recall that So = {ei, . . . , e^} CM.'^ denotes an orthonormal system. 

Lemma 6.2. Assume that = 1^ and that the condition ^/{p, 6, S, 3X) holds with some 
< p < 1 and S > 0. Write m = \pN/{5s)~\. Then, for any 0<A<B,bEW^ and 
'J > 0, we have 

B 



1 1 ^b{t) I ^ < ^ + C7(^) ipNr log 4 ' (6-1) 



with 



A 

B 



sup sup / vW4)^, ^{t)= Ee{tQ[Y + b]} \ (6.2 



A 



where Y = Ui + ■ ■ ■ + Um denote a sum of independent (non-i.i.d.) vectors, and supp is 
taken over all{C{Uj): 1 < J < m} C r((5;D,5). 

Lemma [6.21 is an analogue of Corollary 6.3 from BG (1997a). Its proof is even simpler 
than that in BG (1997a). Therefore it is omitted. 

Lemma 6.3. Assume that = 1^ and that the condition M{p, 5, 5, DX) holds with some 
< p < 1 and (5 > 0. Let 

n= IpN/ilQs)] > 1. (6.3) 
Then, for any 0<A<B,beW'- and 7 > 0, 

B B 

J I Mt) I ^ < c,{s) (pN)'^ log ^+snp J ^/Ee{t{QW,W')/2} , (6.4) 



A 



and for any fixed t G 



^fe(t)| < c^{s){pN)"' + suv^E e{t{QW,W')/ 2], (6.5) 

r ^ 
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where W = Vi + ■ ■ ■ + Vn and W = V( + ■ ■ ■ + are independent sums of independent 
copies of random vectors V and V respectively, and the supremum supp is taken over all 
£(V),£(ner(5;D,5). 

Note that this lemma will be proved for general iS, but in this paper we need S = So 
only. Moreover, a more careful estimation of binomial probabilities could allow us to 
replace c^(s) (pN)~'^ in (16. ip . (16. 4p and (16. 5p by c(s) exp{— cpA^} (see e.g. Nagaev and 
Chebotarev (2005)). However, we do not need to use this improvement. 

Proof of Lemma \ 6.3[ Inequality (16. 5 p is an analogue of the statement of Lemma 7.3 
from BG (1997a). Its proof is even simpler than that in BG (1997a). Therefore it is 
omitted. 

Let us show that 

B B 

1 1 Mt) I ^ < c,{s) {pNy-y log ^+snp J ^JEe{t{QW,W')/2} ^ , (6.6) 

A A 

where W = Vi-\ — ■ + Vn and W = V( + ■ ■ ■ + V^ are independent sums of of independent 
(non-i.i.d.) vectors, and supp is taken over all { £(V^), £(Vj') : 1 < j < C r(5; D, 5). 
Comparing (16. 4p and (16. 6p . we see that inequality (16. 6 p is related to sums of non-i.i.d. 
vectors {V,} and {V^'} while inequality (16. 4p deals with i.i.d. vectors. Nevertheless, we 
shall derive (16. 4p from (16. 6p . 

While proving (16. 6 p we can assume that pN > Cg with a sufficiently large constant Cg, 
since otherwise (16. 6p is obviously valid. 

Let ip{t) be defined in (16. 2p . where Y = Ui + ■ ■ ■ + Um denote a sum of independent 
(non-i.i.d.) vectors with {C{Uj) : 1 < j < m} C r(5;©,>S), m = \pN/{hs)']. 

We shall apply the symmetrization Lemma 16. 1[ Split Y = T + Ti + T2 into sums 
of independent sums of independent summands so that each of the sums T, Ti and T2 
contains n = [pA^/(16s)'| independent summands Uj. Such an n exists since pN > Cs 
with a sufficiently large c^. Lemma [6.11 implies that 

2(^(t) <Ee{2t(Qf,fi)}+Ee{2t(Qf,f2)}. (6.7) 

Inequality (16.61) follows now from (16. 7p and Lemma [6. 2[ 

Let now W = Vi + ■ ■ ■ + Vn and W = V-[ + ■ ■ ■ + be independent sums of of 
independent (non-i.i. d.) vectors with {C{Vj),C{V-) : 1 < j < n] C r{6;B,S). Using 

that all random vectors Vj are symmetrized and have non-negative characteristic functions 
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and applying Holder's inequality, we obtain, for each t, 

n 

Ee{t{QW,W')} = E^,^l[Ey^e{t{QV„W')}) (6.8) 
< [llE^,{Ey^e{t{QV,,W')}y) " (6.9) 

n , , 

= (^l[E^,{E;p^e{t{QT,,W')})) (6.10) 

/ 1 . ~ \ 1/" 

= [l[Ee{t{QT„W')}) , (6.11) 

i=i 

where Tj =^ ^JLi ^ji denotes a sum of i.i.d. copies Vji of Vj which are independent of all 
other random vectors and variables. ^ ^ 

Repeating the steps fl6.8p -( l6?TT]) for each factor E e{ t (Q Tj, W')} instead of the expec- 
tation E e{t {QW, W)} on the right hand side separately, we get (with Tj, =^ J2?=i ^kii 
where V^i are i.i.d. copies of VI independent of all other random vectors) 



2 



. / - - - - ^ ^ \ l/n~ 

Ee{t{QW,W')}<(J{\[Ee{t{QT^X)}) ■ (6-12) 

j=i k=i 

Thus, using f l6.12p and the arithmetic-geometric mean inequality, we have 

B ^ n n 

B 



S ^EE/(Ee{((Qf,.n>/2})"'^ 



j=i k=i 

B 



< sup j ^Ee{t{QT,T>)/2} , (6.13) 



where T = Ui + ■ ■ ■ + Un and T' = U[ + ■ ■ ■ + are independent sums of independent 
copies of random vectors U and U' respectively, and the supremum supp is taken over all 
C{U), C{U') e r((5;©, iS). Inequalities (16.61) and (16.131) imply now the statement of the 
lemma. □ 

The following Lemma 16.41 provides a Poisson summation formula. 
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Lemma 6.4. Let Ke z > 0, a,b & W and E> : W ^ W be a positive definite symmetric 
non- degenerate linear operator. Then 

exp|— zS[m + a] + 27rz (m, 6)} 

2 

= {det(S/7i)y^^^ z-'/^ exp{-27ii {a,b)}J2^^p{- + b] - 2m {a,l)y 

where E>^^ : W ^ W denotes the inverse positive definite operator for E>. 

Proof. See, for example, Flicker (1982), p. 116, or Mumford (1983), p. 189, for- 
mula (5.1); and p. 197, formula (5.9). □ 

Let the conditions of Lemma 16.31 be satisfied. Introduce one-dimensional lattice prob- 
ability distributions Hn = >C(^„) with integer valued ^„ setting 

P{^„ = A;} = gxp {-k^/2n} , for A; G Z. 

It is easy to see that An x 1. Moreover, by Lemma [6.41 

Hn{t) > 0, for all t e R. (6.14) 

Introduce the s-dimensional random vector (n having as coordinates independent copies 
of ^n- Then, for m = (mi, . . . , m^) G Z'*, we have 

q{m) = P{Cn = m} = A'^n-'^"^ exp {-\\mf /2n} . (6.15) 

Lemma 6.5. Let W = Vi + ■ ■ ■ + Vn and W = + ■ ■ ■ + denote independent sums of 
independent copies of random vectors V and V' such that 

V = eizi^ VEsZs, V' = es+iz[^ V e2sz'^, 

with some Zj,Zj G M.'^. Introduce the matrix = {bij{t) : I < i,j < s} with bij{t) = 
t{Qzi,z'j). Then 

E e{t{QW,W')/A} E e{ (B^Cn, C) } + exp {-cn} , for all t G M, 

where ^'^s independent copies of (n OLnd c is a positive absolute constant. 

Proof. Without loss of generality, we shall assume that n > ci, with a sufficiently large 
absolute constant ci. Consider the random vector Y = (ei, . . . , e,) &W with coordinates 
which are symmetrizations of i.i.d. Rademacher random variables. Let R = [Ri, . . . ,Rs) 
and T denote independent sums of n independent copies of y/2. Then we can write 

E e{t{QW,W')/4} =E e{{EtR,T)}, for all t G M. (6.16) 



EXPLICIT RATES OF APPROXIMATION IN THE CLT 29 

Note that the scalar product (-, ■) in E e{ (B^i?, T) } means the scalar product of vectors 
in W. In order to estimate this expectation, we write it in the form 

Ee{{MtR,T)} = EERe{{M,R,T)} 

= p(fn) p(m) e{ (B^m, m) }, (6-17) 

with summing over m = (mi, . . . , m^) ^rn = (mi, . . . , m^) G V and 

s s 

p{m) = P{R = m} = l[p{R,=m,} = (6.18) 

if max |mJ < n and p{m) = otherwise. Clearly, for fixed T = m, 

i<i<s 

En e{ (Bti?, T) } = ^ p(m) e{ (B^ m, m) } > (6.19) 

is a value of the characteristic function of symmetrized random vector MfR. Using Stir- 
ling's formula, it is easy to show that there exist positive absolute constants C2 and C3 
such that 

P{Rj = TTT-j} ^ n'^^"^ exp |— m^/2n} , for \mj\ < C2n, (6.20) 

and 

P{|i?j| > Can} < exp{-C3n} . (6.21) 
Using fleTrjl - dOTj) . we obtain 

Ee{{MtR,T)} ^ g(m) ^ p(m) e{ (Bt m, m) } + exp {— Csn} 

= p("^) g(m) e{ (B^m, m) } + exp {— csn} 

= E Ec„ e{ {MtR, Q } + exp {-cgn} 

= Ee{(Bii?,C„)} +exp{-C3r2}. (6.22) 
Now we repeat our previous arguments, noting that 

E^„e{(Bii?,C„)} = J2 ^^^) 4{^tR,Tn)}>0 (6.23) 

is a value of the non- negative characteristic function of the random vector (see (16.141) ). 
Using again fl6.20p and fl6.2ip . we obtain 

E e{ {MtR, Cn) } <s E e{ (B^C™, Q } + exp {-cgn} . (6.24) 
Relations fl6.16p . (16.220 and fl6.24p imply the statement of the lemma. □ 
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Let US estimate the expectation E e{ (Bt^„, Q'^ } under the conditions of Lemmas 16.31 
and 16.51 assuming that s = c?, D = C^"*^/^, 5 < l/(5s), n > C4,, where C4 is a sufficiently 
large absolute constant, and 

\\C-'^^hj - ejW < 6, \\C-^/^z'j - ejW < S, for 1 < j < s, (6.25) 

with an orthonormal system S = So = { ei, 62, . . . , } involved in the conditions of 
Lemma 1631 We can rewrite E e{ (B^^^, (!^) } as 

Ee{(BtCn,0} = Yl ^(^) Y qMe{{Mtm,m)}. 

Thus, by ( 16151) . 
Denote 

r = V27rhi. (6.26) 

Applying Lemma |6^ with E> = Ig, z = l/2n, a = 0, b = (271)^-*^ B^m and using 
that y4„ X 1, we obtain 

Ee{(BtCn,0} ^ exp{ -271^72 || / + (27r)"X"^f - ||rrif/2^} 

«,mGZ= 

<s r'-' Y exp{ -r^ ||m - tVmf - ||m||V^^}, (6.27) 

m,mGZ'' 

where V : — )■ M'* is the operator with matrix 

V=(27r)"%. (6.28) 

Note that the right-hand side of ( I6.27P may be considered as a theta-series. 

Denote yk = C~^^'^Zk, 1 < k < s. Let Y be the (s x s)-matrix with entries {ej^yk), 
where index j is the number of the row, while k is the number of the column. Then the 

matrix F =^ Y*Y has entries {yj,yk)- Here Y* is the transposed matrix for Y. According 
to fl6.25p . we have 

WVj - GjW < (5, for 1 < j < s, (6.29) 

Let us show that (cf. BG (1997a, proof of Lemma 7.4)) 

||Y||<3/2 and HY^^ < 2. (6.30) 

Since So = { ei, 62, . . . , Cs } is an orthonormal system, inequalities fl6.29p imply that Y = 
Is + A with some matrix A = {ajj} such that |ajj| < 5. Thus, we have ||A|| < ||A||2 < s5, 
where ||A||2 denotes the Hilbert-Schmidt norm of the matrix A. Therefore, the condition 
6<l/[hs) implies ||A|| < 1/2 and inequalities flCTj) . 



EXPLICIT RATES OF APPROXIMATION IN THE CLT 31 

The matrix F is symmetric and positive definite. Its determinant is the product of 
eigenvalues which (by f l6.30p ) are bounded from above and from below by some absolute 
positive constants. Moreover, 

(detY)^ = (det Y*)^ = detF x, 1 x ||F|| x ||Y||. (6.31) 

Define the matrices Y and F, replacing Zj by z'^ in the definition of Y and F. Similarly 
to (I6.3ip . one can show that 

(det Y)^ = (det Y*)^ = detf x, 1 x ||F|| x ||Y||. (6.32) 

Let G and G be the (s x s)-matrices with entries (cj, Q-Zfc) and (cj, z^) respectively. Then, 
clearly, G = QC^/^Y and G = Y. Therefore, 

Ml = G*G = Y*C^/2q(ci/2y. (6.33) 

Moreover, = 1^ implies that |detQ| = 1 and ||Q|| = 1. Using relations (I6.28p and 
flOB -f lOHll . we obtain 

I det V I X, I det Bi | x, det C, (6.34) 

and 

||V|| < pill < ||C|| < (g^35) 



7. Some facts from Number Theory 

In Section[71 we consider some facts of the geometry of numbers (see Davenport (1958) 
or Cassels (1959)). They will help us to estimate the integrals of the right-hand side of 
inequality (16.27^ . 

Let ei, 62, . . . , Crf be linearly independent vectors in M*^. The set 

d 

^ = {SZ'^J-^i -.rij eZ, J = 1,2,..., (7.1) 

3=1 

is called the lattice with basis ci, 62, . . . , e^. The determinant det (A) of a lattice A is the 
modulus of the determinant of the matrix formed from the vectors Ci, 62, ... , 6^: 

det(A) = |det(ei,e2,...,erf)|. (7.2) 

The determinant of a lattice does not depend on the choice of basis. Any lattice A C M"^ 
can be represented as A = AZ'^, where A is a non-degenerate linear operator. Clearly, 
det(A) = |det A|. 

Let mi, . . . , m; G A be linearly independent vectors belonging to a lattice A. Then the 

set 

A' = l^n^m^- : rij G Z, j = 1, 2, . . . , /| (7.3) 
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is an /-dimensional sublattice of the lattice A. Its determinant det(A') is the modulus of 
the determinant of the matrix formed from the coordinates of the vectors mi, m2, . . . , m; 
with respect to an orthonormal basis of the linear span of the vectors mi, m2, . . . , m^. The 

1 /2 

determinant det(A') could be also defined as det((mj, m^), j = 1, . . . , / ) 

Let F : ^ [0,oo] denote a norm on R'^, that is F{ax) = \a\F{x), for a G M, 
and F{x + y) < F{x) + F{y). The successive minima Mi < • ■ ■ < of F with respect 
to a lattice A are defined as follows: Let Mi = inf |F(m) : m 7^ 0, m G A} and 
define Mj as the infimum of A > such that the set |m G A : F{m) < A} contains 
j linearly independent vectors. It is easy to see that these infima are attained, that is 
there exist linearly independent vectors 61, . . . , 6^ G A such that F{hj) = Mj, j = 1, . . . ,d. 
The following Lemma [7.11 is proved by Davenport (1958, Lemma 1), see also Gotze and 
Margulis (2010). 

Lemma 7.1. Let Mi < ■ ■ ■ < M^ be the successive minima of a norm F with respect to 
the lattice . Denote M^+i = 00. Suppose that 1 < j < d and Mj < b < Mj^i, for 
some b > 0. Then 

#{m = (mi,...,mrf) gZ'^: F{m) < b} V {Mi ■ M2 ■ ■ ■ M,)-\ (7.4) 

Representing A = AZ'^, we see that the lattice Z"^ may be replaced in Lemma [7. II by 
any lattice A C M'^. It suffices to apply this lemma to the norm G{m) = F{Am), m G Z'^. 

Lemma 7.2. Let Fj{m), j = 1,2, be some norms in M.'^ and Mi < ■■■ < M^ and 
Ni < ■ ■ ■ < be the successive minima of Fi with respect to a lattice Ai and of F2 
with respect to a lattice A2 respectively. Let C > 0. Assume that M^ C F2{nk), 
k = 1,2, . . . ,d, for some linearly independent vectors ni, n2, . . . , G A2. Then 

Mk:$>dCNk, k = l,...,d. (7.5) 

The proof of this lemma is elementary and therefore omitted. 
Let ||a;||oo = niaxi<j<rf \xj\, for x = {xi, . . . ,Xd) G W^. 

Lemma 7.3. Let A be a lattice in and let < e < 1 . Then 

e-'#H < ^exp{-£||t;f } ^^e-'^'^j^H, (7.6) 

where if =^ | G A : ||f ||oo < 1 } • 

Proof. The lower bound in (17. 6p is almost evident by restricting summation to the set H. 
Introduce for /i = (/ii, . . . , /i^) G Z*^ the sets 



5/= 



1 1 

/^i - Y , /Ui + — ) X ■ ■ ■ X 



1 1 



such that M'^ = U^t -^/i- For any fixed w* e = [w e A^} B^^ we have 

w ~ w* E H, for any w G H^. 
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Hence we conclude for any G 

< i^H. (7.7) 
Since x implies ||x||oo > ||/w||oo/2, we obtain by (17. 7p 

^exp{-£||t;f} < ^exp{-e ||t;||^} 

<d #if ■ ^exp{-£||/i||^/4} 

«, (7.8) 

This conclude the proof of Lemma 17.31 □ 

It is easy to see that Lemma [7.31 implies the following statement. 

Corollary 7.4. Let K he a lattice in and let Cj{d), j = 1,2,3,4, be positive quantities 
depending on d only. Let F{-) be a norm in M.'^ such that F{-) x^^ || ■ ||. Then 

^exp{-ci(rf)||t;f } X, 5^exp{-C2(rf)(F(t;))2} 

veA veA 

-d #{^^ e A : ||f|| < C3(rf)} 

X, #{t;e A : F{v)<c,{d)}. (7.9) 

The proof of Corollary 17.41 is elementary and therefore omitted. Note only that 

i^{v e A : F{v) < X} = i^{v e fi-^A : F{v) < X/jj}, forA,/i>0. (7.10) 
For a lattice A C M*^ and 1 < I < d, we define its archaracteristics by 

a; (A) =^sup{ |det(A')| ^ : A' C A, Z-dimensional sub lattice of A}. (7-11) 

Denote 

a(A) = maxa;(A). (7.12) 

l<l<d 

Lemma 7.5. Let F{-) be a norm in such that F{-) x^^ ||-||. Let c{d) be a positive 
quantity depending on d only. Let Mi < ■ ■ ■ < Md be the successive minima of F with 
respect to a lattice A C M'^. Then 

ai{A) X, (Ml ■ M2 ■ ■ ■ Mi)-\ l = l,...,d. (7.13) 

Moreover, 

«(A)x,#{t;G A : \\v\\ <c{d)}, (7.14) 

provided that Mi <^d 1- 
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For the proof of Lemma 1731 we shall use the following lemma formulated in Proposition 
(p. 517) and Remark (p. 518) in A.K. Lenstra, H.W. Lenstra and Lovasz (1982). 

Lemma 7.6. Let Mi < ■ ■ ■ < Md be the successive minima of the standard Euclidean 
norm with respect to a lattice A C M.'^. Then there exists a basis ci, 62, ■ ■ ■ ,ed of A such 
that 

Mi^dWeil l = l,...,d. (7.15) 

Moreover, 

d 

det(A) J]||e,||. (7.16) 

1=1 

Proof of Lemma \7.5\ According to Lemma [721 we can replace the Euclidean norm || ■ || 
by the norm F{-), in the formulation of Lemma 17.61 Let A' C A be an arbitrary 
/-dimensional sublattice of A and Ni < ■ ■ ■ < Ni he the successive minima of the 
norm -F(-) with respect to A'. It is clear that Mj < Nj, j = 1,2, . . . ,/. On the other 
hand, Mj = F{mj) for some linearly independent vectors mi,m2, . . . ,mi G A. In the 
case, where 

I 

A' = {J^'^jruj : G Z, J = 1, 2, . . . , /}, (7.17) 

i=i 

we have Nj = Mj, j = 1, 2, . . . , Z. In order to justify relation f l7.13p it remains to take into 
account definition (17. lip and to apply Lemma [7]6l Relation fl7.14p is an easy consequence 
of fl7:T3|) . Lemma O and Corollary El □ 



8. From Number Theory to Probability 

In Section m we shall use number-theoretical results of Section [7] to estimate integrals 
of the right-hand side of fl6.27p . Recall that we have assumed the conditions of Lemmas 16.31 
and 16. 5[ s = d, B) = C~^/^, 5 < l/(5s), n > C4, and f l6.25p . for an orthonormal system 
S = So- The notation SL((i, M) is used below for the set of all {d x (i)-matrices with real 
entries and determinant 1. 



Introduce the matrices 



)eSL(2s,M), r>0, (8.1) 



def 



s 

L 



Is -tls 



= ( ^ ) , t G M, (8.2) 

Ut = ( ) e SL(2s,M), t G M, (8.3) 
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and the lattices 

aS(;-_«;)z-. (8.4) 

A,=D,D,-,A= (^^J- ^-"*^Jz^ j = l,2,..., (8.5) 

where 

Vo = af2v (8.6) 

and the matrix V is defined in f l6.28p . Below we shall use the following simplest properties 
of these matrices: 

]D)aBb = Bab, lJaUb = Ua+b and D„U6 = U,2bD„, fora,6>0. (8.7) 

Let Mj^t, j = 1,2, ... ,2s, he the successive minima of the norm || ■ ||oo with respect to 
the lattice 



„ def / ris -rtY \ ^2s 



I?'. (8.8) 



Moreover, simultaneously, Mj^t are the successive minima of the norm F*{-) defined for 
{m,m) G M^'^, m,m E W, by 

F*((m,m)) = max{||m|U, a? \\N~^rri\\J^ (8.9) 
with respect to the lattice 

^t=f'^^^ ^_TiT.,)z2^=D,U„A, where «'=V^. (8.10) 



Using Lemmas 17.21 and [7^ and the equality det(Sf) = 1, it is easy to show that 

Mi,t L (8.11) 

Let M*^ be the successive minima of the Euclidean norm with respect to the lattice Vtt. 
Note that, according to f l6.35p and (18. 9p . 

IHI <,F*(-). (8.12) 

Using fl8.12p and Lemma I7.2[ we obtain 

Ml,<t:sMj,t, 3 = I,..., 2s. (8.13) 

According to Lemma [7. 5[ 

a{Et) a{nt). (8.14) 

Let us estimate a(fit) assuming that r > 1 and (for a moment) t = af^r^^. In this 



case 



= { ) (8.15) 



By relation f l7.14p of Lemma I7.5[ we have 

aiilt) X, G a : < 1/2} = (8.16) 
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where 

/sT = jt; = (m,m) e Z^' : m,m G Z", ||rm - Vomf + ||r~^ Vomf < 1/4 }. (8.17) 

Let us estimate from above the right-hand side of fl8.16p . If f = (m,m) G K, then 

1 r 

r ||m|| < \\rm — Vom|| + ||Vom|| < y + y — (8.18) 

Hence m = and ||Vo?7T|| < 1/2. It remains to estimate the quantity 

=^#{mG Z' : II Vom|| < 1} > (8.19) 

Let A^i < • ■ ■ < be the successive minima of the Euclidean norm with respect to 
the lattice YqIj" . Let ei, 62, . . . , be the standard orthonormal basis of 1/. By f l6.35p 
and (18. 6p . we have ||Voej|| < 1, j = 1, 2, . . . , s. Therefore, using Lemma [7.2[ we see that 
iVi < • ■ ■ < iV, < 1. By f lCTj) . dHU), (I8A9|) and by Lemmas EH O and EH 

R X, (iVi ■ iV2 ■ ■ -iV,)-! X, (det Vo)-^ x, af (detC)-\ (8.20) 

Hence, using dHUD, dHUD and <K7^ we conclude that 

a{VLt) <s (detC)"\ for r > 1 and t = a^^r-^. (8.21) 

Let now t G M be arbitrary. By ([83]), f lHTTTj) . flHTHjl . Lemmas Ell Ol and CorollaryEl 

exp| — ||m — t Vm||^ — ||m||^/r^ } = exp| — ||f |p } 

#{i;GSt : ||i;|| < 1} 

a(Ht) a(a). (8.22) 

Now, by flOTj) . f lHJOj) and (E22]), we have 

Ee{(BtCn,0} =^"'"(DrU„A), where M = ait. (8.23) 

Let us estimate the quantity i?*, t G M, defined in f l8.22p assuming that r > 1 and 
|rt| < c* 0"]"^, where c* > 1 is an arbitrary quantity depending on s only. By Corollary Ell 
we have 

Rt -s #Ko, (8.24) 

where 

Ko = {v = (m,m) G Z^" : m,m G Z^ ||rm-rtVm||2 + Hr'^mf < (20*)"^}. (8.25) 

If w = (m,m) G Ko, r > 1 and |rt| < c* (T^^ then, by (I05l) and (ESHD, 

1 r 

r ||m|| < ||rm — rtYfn\\ + |rt| ||Vm|| ^ ^ + y — (8.26) 
Hence m = and |rt| ||Vm|| < (2c*)~^ < 1. It remains to estimate the quantity 

S = i^{m e : |rt| ||Vm|| < 1} > #7^0- (8.27) 

Let -Pi < ■ ■ ■ < -Ps be the successive minima of the Euclidean norm with respect to the 
lattice |rt| VZ^. Let 61,62, ■■■ ,es be the standard orthonormal basis of Z**. By f l6.35p . 
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we have Vej|| <^s 1, j = l,2,...,s. Therefore, using Lemma 17.21 we see that 
Pi < ■ ■ ■ < Ps -^s I- By (Km . f lOTD and Lemmas O and EH 

5 X, (Pi • P2 ■ ■ • Ps^^ -s (det(|rt| V))-^ x, \rtr (det C)-^ (8.28) 
Hence, using IK2M . f lOTll and IK2M . we conclude that 

Rt <s \rt\~" (detC)-\ for r > 1 and \rt\ < cla^^. (8.29) 
Now, by flOTj) . (E22]) and (ESHD, we have 
Ee{(BtCn,C)} «. r-^^i 

r"^' |tr' (detC)"\ for r > 1 and \rt\ < c*a^^. (8.30) 

It is easy to verify that 

-2 -1 

P ' |tr^ (detC)-i— <,r-Vl^(detC)-^/^ (8.31) 

for any Cs depending on s only. Note that al (detC)"^/^ > 1. Using fl8.23p . fl8.3ip and 
Lemmas 16.31 and 16. 5[ we derive the following lemma. 

Lemma 8.1. Let the conditions of Lemma 16.31 be satisfied with s = d, 1} = C^^^^, 
^ < l/(5s) and with an orthonormal system S = So = {ei, . . . ,6^} C . Let Cg be an 
arbitrary quantity depending on s only. Then, for any b and r > 1, 

I $,(t/2) I ^ {pN)'^ a{ (det C)'^'^ + r'^'^ sup I («(©.U, A))^/' — 



(8.32) 

where r, a(-), D,. Ut and t/ie lattice A are defined in relations (16. 3p . (I6.26p . (I6.28p . (17. lip . 
(I7.12p . (18. ip . (18. 3p and (18. 4p and m Lemma 1^31 T/ie supr means here the supremum over 
all possible values of Zj,Zj G M"^ {involved in the definition of matrices and V) swc/i 
that 

\\C-^^^Zj - ejW < 5, \\C-'^^^z'j - ejW < 6, for 1 < j < s. (8.33) 

Moreover, for any b E W^, r > 1 and 7 > and any fixed t G M satisfying \rt\ < c* ajf^, 
where c* > 1 is an arbitrary quantity depending on s only, we have 

I ^bit) I (pN)-^ + r-' \t\-'/^ (det C)-^/2_ (g_34) 

Let V = (m, m) G M^^, m,fn and t G R. Then 

m + tm = (1 + t^) m + t (m — tm). (8.35) 
Equality (I8.35P implies that 

||m + tm|| ll^^ll + ll'^ — ^?^||, for |t| 1. (8.36) 

Hence, 

r ||m — tm|| + r"""" ||m + tm|| r ||m — tm|| + r"""" ||m||, for r ^ 1, |t| 1- (8.37) 
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According to fl8.1l) - fl8.3l) . we have 

3rVtV = {r{m — tWi), r~^m) and 3rKtV = {r(m — tm), r'^ {fn + tm)). (8.38) 

It is clear that the operators DrU^ and D^Kt are invertible. Therefore, using fl8.37p 
and (18.381) and applying Lemmas 17.21 and I7.5[ we derive the inequality 

a(©,.U^^]) a(D^Ktfi), for r > 1, |t| 1, (8.39) 

which is valid for any lattice Q C M^*. 

Let T be the permutation (2 s x 2s)-matrix which permutes the rows of a (2 s x 2 s)- 
matrix A so that the new order (corresponding to the matrix TA) is: 

l,s + l,2,s + 2, ...,s,2s. 

Note that the operator T is isometric and A i— t- AT^^ rearrange the columns of A in the 
order mentioned above. It is easy to see that 

ajirn) = aj{n), j = l,...2s, and a{Tn) = a{n), (8.40) 

for any lattice f2 C M^^. 
Note now that 

TD^KtAj = TD,,KtT-^TAj = W^A^, (8.41) 
where Aj is a lattice defined by 

Aj = TAj (8.42) 

and where Wt is (2s x 2s)-matrix 







O2 : 


O2 \ 




O2 


Gr,t '■ 


O2 




O2 


O2 : 


^r,t j 



(8.43) 



constructed of (2 x 2)-matrices O2 (with zero entries) and 

g./=(a ^'-''^ 

Let \t\ < 2 and 

6 = arcsin(t (1 + fy^^^) or, equivalently, t = tan6'. (8.45) 

Then we have 

1^1 < c* = arcsin(2/y5), cos0 = (1 + t')-^/^ sinO = t {1 + t^y^^^ (8.46) 
It is easy to see that 

G^^, = (1 + ^2)1/2 p^^^ ^g j^Y) 

and 

Wt = (1 + t2)i/2©,Ke, (8.48) 
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where 



\ O2 O2 

are (2 s x 2 s)-matrices with 
dci f r 



and Kfl 



O2 Ke 



h \ 

h 



\ O2 O2 
dcf / COS 6 — sin ( 



(8.49) 



' _i and Ke =\ . . G SL(2, 

r ^ / V smfc' cosy / ^ 



.50) 



Substituting ( 18.48^ into equahty (I8.4ip . we obtain 

TB)rKtKj = (1 + t^f^ Aj. (8.51) 
Below we shall also use the following crucial lemma of Gotze and Margulis (2010). 

Lemma 8.2. Let Kg and 

/ H O2 : O2 \ 



^.52) 



\ O2 O2 : M / 

be {2dx 2d) -matrices such that M e G = SL(2,]R) and Kg is defined m (I8.49p and f l8.50p . 
Let 13 he a positive number such that (3d > 2. Then, for any EI G ^ and any lattice 
A c M^fiJ^ 

f.27r 

■.Psin / _ / A ||ff||/3d-2_ (8.53) 



(a(eKeA))'^ci^ <^,rf (a(A)) 
Rere ||EI|| zs t/ie standard norm of the linear operator EI : — t- M^. 



Consider, under the conditions of Lemma 18.11 



def 



^-2^_2+4/3/2 



, ^ \ dt 



®*(*/2)|4^. 



c. cr, r 



By Lemma [8. H we have 



where 



with 



Jo <s {pN)-^ al (det C)-i/2 + r"/^ sup J, 



i=2 



dcf 



0-1)-^ 



(a(D.U,A))^/^^, j = 2,3,...,p'=^M+l. 



154) 



^.55) 



(8.56) 



(8.57) 
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Changing variable t = vj""^ and v = w + j in Ij and using the properties of matrices 
and U^, we have 

(«(D,U^,-.U,-iA))'/'^. (8.58) 
By (EB, 

= D^j-i = ro^j-i U^Dj. (8.59) 
According to (ESHD and (l839ll . 

(a(D,j-iUtAj))^^^cit, (8.60) 

where the lattices A^ are defined in (ES]) (see also (JEll), (ES} and (El). Using (1831) . 
(18J[5|1 and (ESID, we see that 

a(Aj) (T?. (jg^q-i_ ^g_g^) 
By ([839]), fl8:iOD and fl83T|) . we have 



^lUfAj) a(D^j-iKfAj-) = «(TDri-iKtAj) 



a(D,j-iKe Aj), (8.62) 

for \t\ <5 1, r > 1, j = 2,3, . . . ,p, where and 9 are defined in and flOHl) 

respectively. Using fISllHD . (EMj), (EIHD, (181621) and Lemma ES (with ci = s), we obtain 

,1/2 r* A ^^l/2 



/ {a{Brj-iVthj)Y''^ dt / (a(D^j-iKe Aj-)) 



cos^ 6 

27r 



/•ZTT 

Jo 



if s > 5. It is clear that = rj'K Therefore, according to IK^ . (K^ . (ESH]) 

and 



-{rr'r/'-'{a{A,)Y^'. (8.64) 



y(rj~')^/'-'(«(^ NNi/2 

By ( ]836|1 . (ESU) and (ESID, we obtain, for s > 5, 

p -I 

J (det C)-^/2 ^ A (rj-i)^/2-2 r^/^-^ ^^^^ ^)-i/2_ ^g g^^ 



i=2 



EXPLICIT RATES OF APPROXIMATION IN THE CLT 



41 



By dnSD, fl835|) and (1835|) . we have r x, (A^p)i/2 and 



Jo <s (det C)-'/^ (iVp)"' ^1 (det Cy^^^. 



(8.66) 



It is clear that in a similar way we can establish that 




$6(t/2) I ^ (det C)-^/2 (jyp)-i (^g^ 



1/2 



(8.67) 



for any quantity c(s) depending on s only. The proof will be easier due to the fact that t 
cannot be small in this integral. 

Thus, we have proved the following lemma. 

Lemma 8.3. Let the conditions of Lemma 16.31 be satisfied with s = d > 5, B) = C^^^^, 
S ^ 1/(5 s) and with an orthonormal system S = So = {ei, . . . ,es} C Mf^. Let ci{s) and 
C2(s) be some quantities depending on s only. Then there exists a Cg such that 
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